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It has been known that an anti-unitary symmetry such as time-reversal or charge conjugation 
is needed to realize Z 2 topological phases in non-interacting systems. Topological insulators and 
superconducting nanowires are representative examples of such Z 2 topological matters. Here we 
report the first-known Z 2 topological phase protected by only unitary symmetries. We show that 
the presence of a nonsymmorphic space group symmetry opens a possibility to realize Z 2 topological 
phases without assuming any anti-unitary symmetry. The Z 2 topological phases are constructed in 
various dimensions, which are closely related to each other by Hamiltonian mapping. In two and 
three dimensions, the Z 2 phases have a surface consistent with the nonsymmorphic space group 
symmetry, and thus they support topological gapless surface states. Remarkably, the surface states 
have a unique energy dispersion with the Mobius twist, which identifies the Z 2 phases experimentally. 
We also provide the relevant structure in the K-theory. 

PACS numbers: 


Introduction .— Symmetry is a key for recent develop¬ 
ments on topological phases. For instance, time-reversal 
symmetry and its resultant Kramers degeneracy are es¬ 
sential for the stability of quantum spin Hall states [1, 2] 
and three-dimensional (3D) topological insulators [3-5]. 
Also, the particle-hole symmetry (or charge conjugation 
symmetry) in superconductors makes it possible to re¬ 
alize topological superconductors [6-18] which support 
exotic Majorana fermions on their boundary. Based on 
these symmetries, many candidate systems for topolog¬ 
ical insulators and superconductors have been proposed 
theoretically and examined experimentally [19-24]. 

In addition to the general symmetries of time-reversal 
and charge-conjugal ion, materials have their own sym¬ 
metry specific to the structures. In particular, crys¬ 
tals are invariant under space group symmetry, like in¬ 
version, reflection, discrete rotation and so on. Such 
crystalline symmetries also provide a new class of topo¬ 
logical phases, which are dubbed topological crystalline 
insulators [25, 26] and topological crystalline supercon¬ 
ductors [27-30]. Surface states protected by crystalline 
symmetry have been confirmed experimentally [31-33]. 
Furthermore, a systematic classification of such topo¬ 
logical phases and topological defects has been done 
theoretically [34-36]. 

In the study of topological crystalline insulators and 
superconductors, much attention has been paid for those 
protected by point group symmetries[37-39]. However, 
point groups are not only allowed crystalline symme¬ 
tries. Space groups contain a transformation which is 
not a simple point group operation but a combination 
of a point group operation and a nonprimitive lattice 


transformation. This class of transformations is called 
nonsymmorphic. In spite that many crystals have such 
nonsymmorphic symmetries, only a few has been known 
for their influence on topological phases [40, 41]. 

In this paper, we show that the presence of nonsym¬ 
morphic space group symmetries provides unique Z 2 
topological phases. Being different from other known Z 2 
phases, the new Z 2 phases need no anti-unitary symme¬ 
try like time-reversal or charge-conjugation. We present 
the Z 2 topological phases in various dimensions, which 
are closely related to each other. In two and three dimen¬ 
sions, the Z 2 phases may have a surface consistent with 
the nonsymmorphic space group symmetry, and thus 
they support topological gapless surface states. Unlike 
helical surface Dirac modes in other Z 2 phase, the sur¬ 
face states have a unique energy dispersion with Mobius 
twist, which provides a distinct experimental signal for 
these phases. The Z 2 topological stability of the surface 
states and a relevant strucuture in the AT-theory are also 
discussed. 

Nonsymmorphic chiral symmetry in ID — As the sim¬ 
plest example, we first consider a ID system. In one- 
dimension, no nonsymmorphic operation is consistent 
with the existence of a boundary, and thus no bound¬ 
ary zero energy state is topologically protected by this 
symmetry. Nevertheless, we can show that an interest¬ 
ing non-trivial Z 2 bulk topological structure appears by 
a nonsymmorphic unitary symmetry. The ID system is 
also useful to construct Z 2 nontrivial topological phases 
in higher dimensions, which have gapless boundary states 
protected by nonsymmorphic symmetries. 

The symmetry we consider is a nonsymmorphic version 
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FIG. 1: (a) Two inequivalent sites A and B in the unit cell. 
(b)Topologically different trajectories {x{kx),y{kx))- (c) In¬ 
sulating state I (top) and II (bottom). 


of the chiral symmetry: In stead of the ordinary chiral 
symmetry, 


{r,iJiD(fc,)} = 0, r^ = i, (1) 


where F is given by a fc 2 ;-independent unitary matrix, we 
consider a fc^-dependent chiral symmetry with 

{riD(fc,),iFiD(fcx)} = 0, . (2) 

By imposing 27r-periodicity in kx on T{kx), the simplest 
^luikx) is 

riD(fc.) = , (3) 


where riD(A:a;) acts on two inequivalent sites A and B 
in the unit cell. As illustrated in Fig.1(a), riD(A:a;) ex¬ 
changes these two sites, followed by a half translation in 
the lattice space. 

The Hamiltonian with the nonsymmorphic chiral sym¬ 
metry has a generic form 


Hiu{kx) 


( x{kx) -iy{kx)e 
\iyikx)e'^’^^/‘^ -x(kx) ) ’ 


( 4 ) 


with real functions x{kx) and yikx)- The 27r-periodicity 
of the Hamiltonian, HiD^kx + 27r) = HiD^kx), implies 


x{kx + 2TT) = x{kx), y{kx+ 2 tt) =-y{kx). (5) 


Because the eigenvalues of the Hamiltonian are E(kx) = 
±\/fs(A(^))pT7y(^T)P, the system is gapped at i? = 0 
unless the vector (a;(fcx), y{kx)) passes through the origin 
(0,0) at some kx- 

Now we will show that the Hamiltonian (4) has two 
distinct topological phases: As we show in Fig. 1(b), the 
Hamiltonian defines a trajectory of {x(kx),y{kx)) in the 
xy-plane, when kx changes from 0 to 27r. From the con¬ 
straint of Eq.(5), the trajectory forms an open arc, not 
a closed circle, and the end point (a;(27r), y(27r)) must be 
the mirror image of the start point (a;(0),y(0)) with re¬ 
spect to the i-axis. The open trajectory passes the cc-axis 


odd number of times. More precisely, we have two differ¬ 
ent ways to across the x-axis; if the trajectory pass the 
positive x-axis odd (even) number of times, then it must 
pass the negative x-axis even (odd) number of times. See 
trajectories To and Te in Fig. 1(c). These two different 
trajectories cannot be continuously deformed into each 
other without gap-closing, since they cannot across the 
origin without gap-closing, as mentioned in the above. 
Therefore, by counting the parity of times the trajectory 
passes the positive x-axis, we can identify the two dis¬ 
tinct phases of the Hamiltonian (4). The Z 2 nature of 
the topological phase is discussed in details in Ref. [42]. 

If the parity is odd (even), then the Hamiltonian is 
adiabatically deformed into the fca;-independent Hamil¬ 
tonian Ho {He) in the below, without gap-closing. 

Ho — , He — ^zi (6) 

with the Pauli matrix ai {i = 0,x,y,z). These Hamilto¬ 
nians suggest a simple physical realization of the nonsym¬ 
morphic chiral symmetry. Consider a periodic potential 
with two different local minima A and B in the unit cell. 
See Fig.1(c). If the energy of the local minimum A (B) is 
much higher than B’s (A’s) and tunnelings between local 
minima are neglected, we have an insulating phase I (H) 
in the half Filing, which effective Hamiltonian is given 
by He {Ho). Our argument above implies that these in¬ 
sulating phases are topologically distinct and they are 
separated by a topological quantum phase transition as 
long as one keeps the symmetry (2). Such a periodic 
system could be artificially created by cold atoms. 

Nonsymmorphic Z 2 symmetry in 2D — Much more in¬ 
teresting Z 2 topological phases protected by nonsymmor¬ 
phic symmetries appear in two and three dimensions. In 
these dimensions, a class of nonsymmorphic symmetries 
are consistent with the presence of a surface, and thus 
the symmetry protected gapless edge states may appear. 
Here we present a 2D 2 ' 2 -topological nonsymmorphic in¬ 
sulator, which supports a unique edge state. 

To obtain the Z 2 phase, we use a Hamiltonian map 
that increases the dimension of the system. This map 
keeps the topological structure by shifting symmetries, 
and is known to be useful to classify the topological (or 
topological crystalline) insulators/superconductors [36, 
43]. In particular, the periodic structure of the topo¬ 
logical table is explained by this map. The details of the 
map in the present case and the relevant structure in the 
K-theory are given in Ref. [42]. 

From the Hamiltonian mapping, we obtain a represen¬ 
tative Hamiltonian of a 2D Z 2 topological nonsymmor¬ 
phic insulator, 

H 2 D{kx,ky) = {m + cos ky)Ty (g) Hiu{kx) + sinkyTx 0 cro,(7) 
which has a fca;-dependent nonsymmorphic symmetry 

[U{kx),H 2 D{kx,ky)] = 0, U{kx) = Tx^Tiu{kx), (8) 











and the additional chiral symmetry, 
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{r,i?2D(fcx,fcy)} = 0, r=Tj;0(7o (9) 

where Ti {i = 0,1, 2,3) is the Pauli matrix for the de¬ 
grees of freedom on which P acts. These two symmetry 
operators anticommute 

{r,t/2D(fcx)} = o. (10) 

Here note that the nonsymmorphic symmetry U 2 D{kx) 
commutes with H 2 D{kx,ky), although it is constructed 
from riD(fcx) anticommuting with HiD(kx). Whereas 
any terms consistent with the symmetries (8) and (9) 
can be added to the Hamiltonian, the basic topologi¬ 
cal properties can be captured by Eq.(7). For a gapped 
Hm{kx), the system has a gap unless to = ±1. Using 
the symmetries (8) and (9), we can define a Z 2 invari¬ 
ant, which is nontrivial (trivial) if —1 < to < 1 (to > 1 
or TO < —1)[42]. Without loss of generality, we assume 
in the following that the parity of Hi^^kx) is even, so it 
is topologically equivalent to = az- 

If we consider a boundary parallel to the x-axis, we can 
keep the symmetries (8) and (9). This boundary supports 
gapless edge states when the system is topological (—1 < 
TO < 1): To demonstrate this, consider a semi-infinite 
system ijj > 0) with the edge at y = 0. Since i7iD(fcx) 
is topologically equivalent to CTz, we first consider the 
spacial case of the Hamiltonian (7) with HiY){kx) = cr^. 
In this particular case, the Hamiltonian H 2 Y}{kx, ky) does 
not depend on kx, and thus the topological edge state 
should be a fc 2 ,-independent zero energy state. The edge 
state can be obtained analytically when the system is 
close to the topological phase transition at to = ±1. Near 
the topological phase transition, say at to = — I, the bulk 
gap is nearly closed at ky = 0, so the low energy physics is 
well-described by the effective Hamiltonian obtained by 
the expansion of Eq.(7) around ky = 0. Then, replacing 
ky with —idy, we have the equation for the edge state 

[{m + l + dyl2)Ty 0 cr^ - idyTx 0 cto] ip{y) = 0, (II) 

with the boundary condition tjj(0) = 0 and ijj{oo) = 0. If 
the system is in the topological side near the transition 
i.e. 6m = to -I- 1 > 0, the equation has two independent 
solutions localized at y = 0 

|iAi) = ® (1) sinh (y-2Sm + ly^ 

IV' 2 ) = 0 e~y sinh (^V-2.5m + Ij/^ (12) 

On the other hand, in the non-topological side {Sm < 0), 
the solutions are diverge, and the edge states disappear. 
A similar result is found near another transition point at 
TO = 1. We have also confirmed numerically the existence 
of the zero energy edge mode for the whole region of 
— 1 < TO < 1. 



FIG. 2: (color online). Schematic illustration of edge states 
with Mobius twist, (a) The red (blue) line is an edge state 
in the eigensector of U{kx) with the eigenvalue u = 

{u = (b) An exchange process of the eigensectors. 

For a general fca;-dependent HlYy{kx)^ the zero energy 
edge states have a fc^-dependent energy dispersion. By 
diagonalizing the mixing matrix {ipi\[5m -I- dy/2)Ty 0 
(Hiuikx) — az)\4’j)j the energy is evaluated as E(ky) oc 
±y(kx)- Then, from the constraint (5), there must be an 
odd number of zeros for y{kx) in kx € [0, 27r], and thus the 
energy dispersion becomes helical E{ky) ~ Ec{kx — fco) 
around each zero /cq, as illustrated in Fig.2 (a). 

Since the Hamiltonian H 2 D{kx,ky) commutes with 
U{kx), the helical dispersion is decomposed into chiral 
and anti-chiral ones, each of which is an eigenstate of 
U{kx). These two chiral dispersions are mapped to each 
other by the chiral symmetry T, because T maps a gap¬ 
less state to another one, reversing the slope of the dis¬ 
persion. Furthermore, they belong to different eigensec¬ 
tors of U{kx), because T exchanges the eigenvalues of 
U{kx) due to {r,U{kx)} = 0. Therefore, these two chi¬ 
ral dispersions stay gapless without mixing, as far as the 
symmetries (8) and (9) are retained. 

Whereas the above edge state has a similarity to helical 
edge modes in quantum spin Hall states, their overall 
structure in the momentum space is completely different: 
As is seen in Fig.2 (a), the present edge state has a unique 
energy dispersion with the Mobius twist, which is never 
seen in other Z 2 phases. This twist occurs due to the 
multivalueness of the eigenvalues u = of U{kx): 

When one goes round in the fca;-direction as kx —^ kx+2TT, 
u changes the sign, so a chiral dispersion in an eigensector 
of U{kx) turns smoothly into to an anti-chiral one in 
another eigensector. 

Another remarkable feature of our edge state is that 
the constituent chiral dispersions can exchange their 
eigensectors of U{kx), as illustrated in Fig.2 (b). This 
means that any pair of helical dispersions is topologically 
unstable: When a pair of helical dispersions exits, we can 
always realize the situation where a chiral dispersion co¬ 
exists with an anti-chiral one in the same eigensector of 
U{kx), by exchanging the eigensectors properly. Thus, 
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we can open a gap of helical dispersion by mixing be¬ 
tween the chiral and anti-chiral ones. 

The arguments above clearly indicate that helical edge 
states in this system has a Z-i stability like helical edge 
states in quantum spin Hall systems, although no time- 
reversal symmetry is required and the mechanism of the 
stability is completely different. 

Glide reflection symmetry in 3D — Finally, we consider 
the system with glide reflection symmetry, 

G{kx)H3D{kx,ky,kz)G~^{kx) = H3u{kx,ky,-kz), 

G\kx) = e-^^% (13) 

The glide reflection G{kx) is the combination of reflec¬ 
tion with respect to the xy-plane and translation along 
the x-axis by a half of the lattice spacing. Since G‘^{kx) 
results in a translation by a unit lattice spacing in the 
x-direction, it provides the non-trivial factor. The 

Z 2 invariant defined by the glide reflection symmetry is 
given in Ref. [42]. 

A representative Hamiltonian with glide reflection 
symmetry is given by [42] 

H 3 D{k) = {m + cos kz + cos ky)Ty iSi Hiuikx) 

+ sin kyTx 0 CTo -I- sin kzTz 0 cto, 

G{kx) = Tx 0Tij:,ikx). (14) 

The 3D system is gapped unless m = ±2, 0. The Z 2 
invariant is non-trivial (trivial) when —2 < m < 0 or 
0 < TO < 2 (to < —2 or TO > 2) [42]. 

A surface perpendicular to the y-axis retains the glide 
reflection symmetry, so it may support a gapless surface 
state protected by this symmetry. For instance, consider 
a semi-infinite 3D system (y > 0) with a surface at y = 0, 
which preserves the glide reflection symmetry. In a man¬ 
ner similar to the 2D system, for the special but topolog¬ 
ically equivalent case with iFiD(^x) = o'z, we can obtain 
the surface state analytically near the topological phase 
transition at to = ±2: For to ~ —2, H 3 Y){k) is well ap¬ 
proximated by 

H 3 B = {m + coskz + 1 - dy/ 2 )Ty 0 tJz 

- 0 Cto -I- sin 0 Cto . (15) 

We find that l"^!) and IV' 2 ) in Eq.(12) with 6 m = m + 
cos kz + 1 satisfy the Schrodinger equation, 

Hso'fiiv) = Ei{kz)'fi{y), (16) 

with Ei(kz) = sinfcz and E 2 {kz) = — sinfcz, respectively. 
When the system is in the topological side near the tran¬ 
sition, i.e. —2 < TO < 0, 6m is positive (negative) at 
kz = 0 {kz = tt). Thus, they meet the boundary con¬ 
dition ipi{0) = 0 and ipi{oo) = 0 near kz = 0, while 
they diverge near kz = tt. This means that they form 
surface states with the linear dispersion E{kz) = Ekz 
near kz = 0, which merge into bulk states near kz, = tt. 


On the other hand, in the topologically trivial side, i.e. 
TO < — 2, 6m is always negative, so {ipi) and IV' 2 ) are no 
longer physical states anymore. A similar analysis works 
for 0 < TO < 2, although the surface states appear near 
kz = TT in this case. 

For a general E[iY){kx), the surface states have a dis¬ 
persion in the fcx-direction, as well as in the fcz-direction: 
Like the 2D case, the two surface modes, |V'i) and |'02), 
are mixed. The spectrum of the surface states becomes 

E{kx,kz) = ± cy’^(kx) + sin^ kz. (c is a constant.) 
From the constraint (5), y{kx) has an odd number of 
zeros, and thus the surface states have the corresponding 
odd number of Dirac cones in the spectrum, as illustrated 
in Fig.3. 

In the glide invariant plane at fcz = A (A = 0, tt) in 
the Brillouin zone, the Dirac cone has helical dispersions 
E ~ ±c{kx — ko) in the fca;-direction. Since E[ 3 D{k) com¬ 
mutes with G{kx) at kz = A, the helical dispersion can 
be divided into two eigensectors of G{kx), which have 
chiral dispersion and anti-chiral dispersion, respectively. 
These two chiral dispersions cannot mix, so a single Dirac 
cone is topologically stable. On the other hand, a pair 
of Dirac cone is topologically unstable: From a process 
similar to Fig.2 (b), the eigensectors can exchange with¬ 
out gap-closing. Therefore, from a similar argument in 
the 2D case, helical dispersions for a pair of Dirac cones 
can be gapped. 

As in the 2D case, the obtained surface state has a 
very unique feature: In the fca;-direction, which is the di¬ 
rection of the translation for the glide, the surface state 
has an energy dispersion with the Mobius twist. Further¬ 
more, along the same direction, the surface state is de¬ 
tached from the bulk spectrum. Indeed, by adiabatically 
changing iFiD(fcx) as FIid = <Jz, the surface state be¬ 
comes completely flat at if = 0 in the fca;-direction. This 
feature is never seen in surface Dirac modes in other Z 2 
phases. Any stable Dirac mode in other Z 2 phases bridge 
the bulk conduction and valence bands in any direction 
in the surface Brillouin zone. This remarkable feature in 
the spectrum can be detected by angle-resolved photoe¬ 
mission, which provides a distinct evidence of this novel 
Z 2 phase. 

Summary — We have revealed that nonsymmorphic 
crystalline symmetries such as glide reflection symmetry 
provide a class of novel Z 2 phases. They are related to 
each other by Hamiltonian mapping, which is justified by 
the K-theory [42]. These Z 2 phases predict remarkable 
surface states that have the Mobius twist in the spec¬ 
trum, which can be detectable experimentally. 

Note added — After this work was finalized, there 
appeared a complementary and independent work [44] 
which has some overlap with our results. 
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FIG. 3: (color online). A surface state protected by glide re¬ 
flection symmetry. The spectrum at = A has a Mobius 
twist in the fc^j-direction: Along the fcj^-direction, the red 
branch with the eigenvalue of G{kx) turns into the blue 

one with the eigenvalue — 
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Supplementary Material 


Hamiltonian mapping 

Here we introduce a Hamiltonian mapping which re¬ 
lates topological phase in different dimensions. A similar 
map has been used in the classification of topological in¬ 
sulators and superconductors defined on a sphere k G S‘^ 
in the momentum space[36, 43]. We generalize the idea 
to insulators with nonsymmorphic symmetries. 

First we review the Hamiltonian mapping used in topo¬ 
logical insulators and superconductors. The map is given 
as follows: If a Hamiltonian H{k) on a d-dimensional 
sphere k G S‘^ has chiral symmetry, {T,H{k)} = 0, with 
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the chiral operator F, then the map is 

H{k,0) = siii0H{k) + cos or, 0G[O,7r], (SI) 
and if not, it is 

H{k, 0) = sin^Ty ® H{k) -\- cos0Tx 0 1 , 0 € [0, tt], (S2) 

where 1 is the unit matrix with the same dimension as 
H{k). Since the mapped Hamiltonian H{k,0) is inde¬ 
pendent of fc G at 0 = 0 and tt, the base space 
{k,0) G X [0, tt] oiHik, 0) can be regarded as a (d-l-1)- 
dimensional sphere 5''^+^ by shrinking to a point at 
0 = 0 and TT, respectively. Thus the mapped Hamilto¬ 
nian H{k,0) is defined on Furthermore, it can be 

shown that the map is isomorphic and thus the original 
H{k) and the mapped H{k,0) have the same topolog¬ 
ical structures. This map relates topological insulators 
in different dimensions, and it enables us to study their 
topological phases systematically. 

Using the above isomorphic map, we can construct the 
2D insulator 

H 2 D{kx,ky) = {m + cosky)Ty (g) Hiuikx) 

+ sin kyTx 0 ao, (S3) 

with symmetries 

[U{kx), H 20 {kxi ky)\ — 0; U{kx) = gr ^ Tid {kx\ 
{T,H2D{kx,ky)} = 0, r = T2(g)CTo, (S4) 

which is topologically nontrivial (trivial) for — 1 < m < 1 
(m > 1 or m < — 1). 

The basic idea is as follows: For Hix^{k) on kx G 
, consider the following two Hamiltonians defined on 
{kx,0) € X [0,7r], 

Hn{kx, 0) = sin 0Ty ® Hij^{kx) + cos Otx (g) (Tq, 
Hi^{kx,0) = sin 0Ty (g) [-iFiD(fcx)] + cosOtx (g) cr^^S) 

which are obtained by the isomorphic map (S2). They 
have the symmetry, 

[Uikx),Hn,Likx,0)]=0, {r,iFR,L(fcx,0)} = O, (S6) 

with U{kx) = Tx ® FiD(fcx) and F = Tz g) (Tq- Since 
Hroikx) and —Hi^(kx) have different Z 2 numbers, ei¬ 
ther Hji{kx, 0) or Hi^{kx, 0), but not both is topologically 
nontrivial. These two Hamiltonians coincide at 0 = 0 
and TT, respectively. Thus, by sewing these two Hamil¬ 
tonians at 0 = 0 and tt, as illustrated in Fig. 4, we can 
obtain a system defined on a two-dimensional torus T^. 
The resultant system has a non-trivial Z 2 number, which 
is obtained as the total Z 2 numbers of H^(kx,0) and 
Hhikx,0)- 

To obtain an explicit Hamiltonian of the system on T^, 
we change the variable 0 as 0 = 7r/2 — fcy in H^{kx,0) 


and 0 = fcy — 7r/2 in Hi^{kx, 0), respectively. For the new 
variable, iJR and iFn have the same form as H 2 T>{kx, ky) 

H2Ti{kx, ky) = coskyTy g Hiuikx) + sinkyTx g tq, (S7) 

where FFl and FFr, are smoothly connected at fcy = 7r/2 
and ky = —7r/2, respectively. Equation (S7) is the Hamil¬ 
tonian of the sewn system. Note that we may adia- 
batically add a term preserving the symmetries (S6) to 
the Hamiltonian without changing its topological prop¬ 
erty unless the bulk gap of the system closes. Thus 
we can finally modify (S7) in the form of Eq.(S3) with 
— 1 < m < 1. 

In a similar manner, we can obtain a system on 
with trivial Z 2 topology. In this case, we use the same 
Hamiltonian for HYi{kx,0) and Hi^(kx,0), 

Hn{kx,0) = Hi^{kx,0) 

= sin0Ty ® Hiiyikx) + cosOtx ® (To, (S8) 

with 0 G [0,7r]. Even when FIr and Hi^ have non-trivial 
Z 2 numbers, they are canceled by sewing them at 0 = 0 
and 0 = TT. An explicit form of the sewn Hamiltonian 
is obtained as follows. Because sin0 > 0, we can add 
a positive constant m to sin0 in Eq.(S8) without gap¬ 
closing, 

Hn{kx,0) = Hi,{kx,0) 

= {m + sin 0)Ty g Hiuikx) + cos0Tx g cto, (S9) 

where we gradually increase m as it satisfies m > 1. 
Then we can adiabatically change the coefficient of sin 0 
in Hi,{kx,0) as sin0 —>• — sin0, without gap-closing. As 
a result, FFr and FIl can be 

Hn.{kx,0) = {m + sin0)Ty g FFiD(fcx) + cos0Tx g ctq, 
Hh{kx,0) = (m- sin0)Ty g HiD{kx) + cos0Tx g o((gl0) 

with m > 1. Finally, by changing the variable 0 as 0 = 
7r/2 — ky in H^{kx,0) and 0 = ky — 7r/2 in Ho,{kx,0), 
respectively, we find that FFr and FFr have the form of 
Eq.(S3) with m > I, where FFr and FFr are smoothly 
sewn up at ky = ±7r/2. We note that if we take the 
starting Hamiltonians as 

FFR(fc„0) =FFR(fc„0) 

= sin0Ty g [-FFiD(fca;)] + cos0Tx g ctq, (Sll) 

we can obtain Eq.(S3) with m < — I, in a similar manner. 
The same idea is available to obtain the 3D insulators 

Hsoik) = {m + cos kz + cos ky)Ty iS) Hiuikx) 

+ sin kyTx g To -I- sin k^Tz g tq, (S12) 

with the glide reflection symmetry, 

G{kx)Hoo{kx,ky,kz)G~^{kx) = Hsuikx, ky,-kz), 

G{kx) = Tx^Tiuikx), (S13) 
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where X{kx) and Y{kx) are N x N hermitian matrices. 
Since Hi-a(kx) is 27r-periodic in k^, Xik^) and Y{kx) 
satisfy 

X{k, + 2tt) = X(fc,), + 27r) = -Y{k,). (S18) 

Now we introduce the following N x N matrix Z{kx) 
Z{k^)=X{k^)+iY{k^), (S19) 

which has the constraint 


FIG. 4: Hamiltonian mapping. Two Hamiltonians 77R(fe,0) Z(fca; + 27r) Z^(kx)- (S20) 

and Hnik, 6) defined on T'^ x fO, ttI are sewn at 0 = 0 and tt. i i . 

Because one can prove the relation 


which is .Z 2 -non-trivial (.Z 2 -trivial) for — 2 < m < 0 or 
0 < m < 2 (m < —2 or m > 2): Since H 21 ) in Eq.(S3) 
is chiral symmetric, we use the isomorphic map (SI) to 
have HB.ikx,ky,6) and Hi,{kx,ky,9), 

HB.{kx,ky,9) = sm9H2j3(kx,ky) +cos9T, 
Hi^{kx,ky,9) = sm9H2j){kx,ky) + cos9r, (S14) 


deti7iD(fca:) = |det2'(fcx)P, (S21) 

detZ{kx) ^ 0 when i7iD(fca:) is gapped at £1 = 0 (namely, 
when deti£iD(fca:) ^ 0). 

Denoting the real and imaginary parts of (leiZ(kx) as 
x{kx) and y{kx), respectively, the relation (S21) implies 
(kx) + y'^{kx) 0 for a gapped i7iD(fcx)- Furthermore, 
from Eq.(S20), we have 


where we denote H 2 D in iJn as different 

from H 2 B in and iJu have the same Z 2 topolog¬ 
ical number as H 2 B and respectively. By jointing 

and Hi^ at 0 = 0 and tt, we can have a Hamiltonian 
Hsd defined on a 3D torus T^. If either or iJn, but 
not both is Z^ 2 -nontrivial, H^d is 2l2-nontrivial. In other 
cases, H^d is .Z 2 -trivial. Then, one can show that with 
a suitable adiabatic deformation, H^d takes the form of 
Eq.(S12) without gap-closing. 

Z2 invariants for nonsymmorphic systems 
ID case 


x{kx+2Tr) = x{kx), y{kx + 2 tt) =-y{kx) ■ (S22) 

Since x(kx) and y(kx) defined here have the same prop¬ 
erty as those in the main text, we can define the Z 2 
invariant in the same manner. 

As is shown in the main text, the simplest Hamil¬ 
tonian with the non-trivial Z 2 invariant is 
which gives {x{kx),y{kx)) = (—1,0). To confirm the Z 2 
nature, consider the direct sum Ho © Ho- In the ba¬ 
sis where riD(fca;) takes the form of Eq.(SI6), Ho © Ho 
gives X{kx) = —l 2 x 2 and Y{kx) = 0. Thus, we find 
{x{kx),y{kx)) = (1,0) for Ho © Ho, which implies that 
Ho © Ho is Z 2 -trivial. 


Here we generalize the Z 2 invariant defined for the 
simplest 2x2 Hamiltonian (4) in the main text, to that 
for the general Hamiltonian. 

The nonsymmorphic chiral symmetry is given by 

{TiD{kx),Hij,{kx)} = 0, r?D(fc,) = e-*'=-. (SIS) 

By imposing 27r-periodicity in k^ on riD(fc 3 ,), a general 
form of riD(£a;) is given by 

/n p-»fcx\ 

riD(fcx) = K Q j © l-AfxAr, (S 16 ) 

with the N X N unit matrix Inxn- In this basis, the 
Hamiltonian HiD{kx) with the nonsymmorphic chiral 
symmetry takes the form 


2 D case 

In this section, we define the Z 2 invariant for the 2D 
Hamiltonian which has the nonsymmorphic symmetry 

[U{kx),H2Ty{kx,ky)] = 0, U{kxf = , (S23) 

as well as the ordinary chiral symmetry, 

{H 2 Dikx,ky),r} = 0, r^ = l. (S24) 

These symmetries are anticommute, 

{r,U{kx)} = 0. (S25) 

Consider the Schrodinger equation given by 

H2D {kx ; ky) \ Un{kx , )) = Efi (^kx , ky^\Un (^kx t ky)) ,(S26) 


^iDikx^ 


( X{kx) -tY{kx)e-^'^^/^\ 
bD(fc,)e*'=-/2 -X{kx) ) 


(S17) 


where n is the band index. We assume that the system 
is gapped at E = 0, and the Fermi energy is inside the 





gap. It is convenient here to use a positive (negative) n 
to represent a positive (negative) energy band. 

Since H 2 D{k^,ky) commutes with U[kx), the solution 
\un{kxTky)) are taken as eigenstates of U{kx) 


that Af{kx,ky) {kx,ky) is a total derivative of a 
function, which yields 


gj[7+(fcx)+j7 (kx)] _ 2 


(S36) 


U{kx)\ut{kx,ky)) = ±e-^’^^^^\ut{kx,ky)). (S27) 


Combining this with Eq.(S34), we finally have 


The chiral symmetry implies that if \u^{kx, ky)) is a pos¬ 
itive (negative) energy band, T\u^{kx,ky)) is a negative 
(positive) energy band. From the anticommutation rela¬ 
tion (S25), it is also found that T\u^{kx, ky)) is an eigen¬ 
state of U{kx) with the eigenvalue Therefore, 

we can place the relation 


\ut{kx,ky)) = T\ut^{kx,ky)). (S28) 

A key character of the nonsymmorphic symmetry 
U{kx) is that its eigenvalues do not have the 

same periodicity as U{kx) itself: They change their sign 
when kx ^ kx + 27r. As a result, \u^{kx,ky)) and 
\u~ {kx + 27r, ky)) have the same eigenvalue of U{kx), sat¬ 
isfying the same Schrodinger equation. Thus, they are 
the same state up to a 17(1) gauge factor, 

\utikx + 2TT,ky))=e^^-^'^^’'^y^u-ikx,ky)). (S29) 

This relation gives a non-trivial relation in Berry phases: 
Introducing the gauge field in the momentum space, 

Af{kx,ky) = i'^{u^{kx,ky)\dkiutikx,ky)), (S30) 

n<0 

we define the Berry phases 7^(fcx) as 

= exp ^ dkyA^{kx,ky)'^ . (S31) 

Since Eq.(S29) implies 

Afikx + 27r, ky) = Af{kx,ky) - ^ dkiOn{kx, fcy)(S32) 

n<0 


the Berry phases satisfy 


(kx+'iTr) _ g-il*{kx) 


(S37) 


Using this relation, we can define the Z 2 invariant in 
the same manner as the ID case: Denoting the real and 
imaginary parts of as x{kx) and y{kx), respec¬ 

tively, we can introduce a nonzero two-dimensional vec¬ 
tor {x{kx),y{kx))- Then Eq.(S37) gives the constraint 


x{kx+2Tr) = x{kx), y{kx + 2 tt) = -y{kx), (S38) 


which is exactly the same as Eq.(5). Therefore, if the 
trajectories {x{kx),y{kx)) passes the positive a:-axis odd 
(even) number of times, the system is topologically non¬ 
trivial (trivial). 

The Z 2 invariant of the Hamiltonian (7) is evaluated 
as follows. It is sufficient to consider the case with 
Hinikx) = (Xz since HiYtikx) can deform into Uz without 
gap-closing. H 2 T){kx,ky) is block diagonal in the diago¬ 
nal basis of U{kx), and in the sector with the eigenvalue 
u = ^e-'’kxl 2 pf U{kx), it is given by 


From this, we obtain 


(S39) 


l^ikx) 


0, for m < — 1 

TT, for — 1 < m < 1 , 

0, for m > 1 


(S40) 


which implies the Z 2 invariant is non-trivial (trivial) if 
— l<m<l (m>lorm<—1). 


3D case 


^i-y+{k^+2Tv) _ gj7 (kx) f dky J 2 n <0 ,ky) ^ (S33) 

From the periodicity in fcy, the integral § dkydkydn 
should be 27 :with a integer and thus we have 

gi7+(fex-l-27i-) _ gi7“(fcx)^ (S34) 

Now we use Eq.(S28). This equation implies 

Aj {kxiky) -\- A^ ikxiky) 

= i'^{ut{kx,ky)\dkiut{kx,ky)), (S35) 

n 

where the summation in the right hand side is taken for 
all n. Therefore, from the completeness relation, we find 


Finally, we define the Z 2 topological invariant associ¬ 
ated with glide symmetry 

G{kx)H3D{kx,ky,kz)G~^{kx) = Hsuikx, ky,-kz), 

G^{kx) = e-^’^-. (S41) 

From solutions of the Schrodinger equation 

H3J){k)\Un{k)) = En{k)\un{k)), (S42) 

we introduce the gauge field Ai{k) in the momentum 

space, 

Ai{k)=i ^ {unik)\dkiUn{k)), (S43) 

E„{k)<Ep 
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FIG. 5: Upper half Brillouin zone. 


where E-p is the Fermi energy. On the glide invariant 
plane at fej = A (A = 0 , 7 r), the glide operator G{kx) 
commutes with H^d, 

[G{kx),H3Bikx,ky,A)]=0, (S44) 

and thus the solutions of the Schrodinger equation are 
simultaneously eigenstates of G{kx), 

G{kx)\ut{kx,ky,A))=±e-^^-/^\u^{kx,ky,A)).{SA5) 

Correspondingly, we can decompose Ai{k) into two parts, 

Ai{kx,ky,T) = Af{kx,ky,A) + A~{kx,ky,A) (S46) 

with 

Af{kx,ky,A) 

= * X! {ut{kx,ky,A)\dk,ut{kx,ky,A)).{S^l) 

En<iE'p 

In a manner similar to U{kx) in 2D, the eigenvalues of 
G(kx) do not have the same periodicity in kx as G{kx) 
itself, and they change their sign when kx ^ kx + 2 tt. As 
a result, we have a twisted boundary condition, 

\utikx + 27T, ky, A)) = (fcx, ky, A))(S48) 

where 9^{kx,ky,A) is a 17(1) phase. 

Now we consider the upper half region of the Brillouin 
zone in Fig.5. From the twisted boundary condition, the 
Berry phases 7 ^(^) along i = a,b,c,d in Fig.5, 

j^^dkyA^{kx,ky,A) (S49) 

satisfy 

7^ (a) =7“(c) (mod. 27r), 

'y~^ {b) ='y~ {d) (mod. 27r). (S50) 

The Stokes’s theorem also leads to 

7(0 + 6)= / Fyzdkydkz (mod. 27r), 

J A 

7^(c—a)= / F^dkxdky (mod. 27r), 

JB 

7 ^( 6 —d)= j F^dkxdky (mod. 27r), (S51) 

Jc 


with 7 (£) = 7 +(£) +7 {(.), Fyz = dkyAz - dk^Az, and 
F^ = dk^A^ — dkyA"^. The modular equality in the 
above equations comes from the ambiguity of the Berry 
phases. 

Using these relations, we find that the following 1 / de¬ 
fines the Z 2 invariant (—l)'^: 

1 

1 / = — 

2 t 

— — 7 "'"(a + 6 ) (mod. 2). (S52) 

IT 

Here note that the modulo-2 ambiguity from the Berry 
phase 7 “'"(a+ 6 ) does not affect on the Z 2 invariant (— 1 )“^. 
In order for (—1)“^ to define the Z 2 invariant, v must be 
an integer. From Eqs.(S51) and (S50), we find that 


Fyzdkydkz + 


IB-C 


Fxydkxdky 


1 

27r 


>B-C 


^xy ^^x dky 


= ^7 (c - a - 6 + d) 

= ^ [ 7 “(c + d) - 7 “(a + 6 )] 

= ^ [ 7 +(a + 6 ) - 7 “(a + 6 )] (mod. 1). (S53) 
Therefore, v is recast into 


V 


1 

27r 



Fyzdkydkz 


7(0 + 6 ) 


(mod. 1)(S54) 


which takes an integer. 

From the formula (S52), we can calculate the Z 2 invari¬ 
ant (—1)*^ for H 3 p,{k) in Eq.(14). Since the Z 2 invariant 
takes the same value unless the gap of the system closes, 
we can choose the special case of Hip,{kx) = cFz- In this 
case, the first and the second terms of the right hand 
side of Eq.(S52) vanish, and thus we only need to eval¬ 
uate 7 ~'"(a + 6 ). On the glide invariant plane at kz = A, 
Hsuikx-, ky, A) is decomposed into in the sector with 
the eigenvalue of G{kx), 

ij^x ; ; A) 

sinfcy z(m + cos A + cosfcy) 

—i(m + cos A + cos/cj,) — sin ky 

(S55) 



From this, we find that 


r 0 , 

for TO < —2 

7+(a) = TT, 

for —2 < TO < 0 

lo> 

for TO > 0 

r 0 , 

for TO < 0 

7+(6) = < TT, 

for 0 < TO < 2 , 

lo> 

for TO > 2 


which implies 

{ 0 , for m < —2 
1 , for —2 < TO < 0 
1, for 0 < TO < 2 ’ 

0 , for TO > 2 


(S56) 


(S57) 
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modulo 2. Therefore, H 3 D{k) in Eq.(14) is topologically 
non-trivial (trivial) if —2 <m< 0 or 0 <m <2 (m > 1 
or m < — 1 ). 


K-theory analysis 

We summarize relevant results in the /tT-theory. Con¬ 
sider a class of nonsymmorphic symmetries, {C/|rx}, 
which consist of a point group operation U accompa¬ 
nying a half translation Tx of the lattice spacing in the 
^-direction. We assume that the point group operation 
[7 is a Z 2 -transformation (namely order-two). The non¬ 
symmorphic symmetry {U\tx} acts on the Bloch Hamil¬ 
tonian H{k) as a fca;-dependent unitary transformation 
Uikx) with U'^ikx) = 

Let us denote the K-group for d-dimensional insu¬ 
lators with the nonsymmorphic symmetry U{kx) as 
Here the superscript {s,t,Tx) identifies the 
symmetries of the insulators: s = 0 , 1 (mod. 2 ) indicates 
the absence (s = 0 ) or the presence (s = 1) of the addi¬ 
tional chiral symmetry. Then t = 0,1 (mod.2) determines 
how the point group operation of Uikx) acts; for s = 0 , 
t specifies the action of U{kx) as [36], 

= { %%% I ;(S58) 

and for s = 1 , 

{Hikx, k),r} = 0, U{kx)H{kx, k)U{kx)-^ = Hikx,k), 

(S59) 


rTJ(h ) — / ^(^a;)r {t — 0 ), 

^ " 1 -t/(fc.)r (t = i). 


where k k represents a Z 2 transformation for 
Finally, Tx represents the half translation in the x- 
direction, as mentioned above. {t,Tx) is an example of 
a twisting of the twisted equivariant iL-theory for topo¬ 
logical insulators and superconductors. [45] 

From the Gysin exact sequence [46] in the A"-theory 
(the twisted version follows from the Thom isomorphism 
theorem [47]), for except in the fca;-direction, the fol¬ 
lowing isomorphism can be shown, 

X S^) 

if U{kx) for x S^) acts on as a global 

symmetry, or 

Ki^^t,r^){Td X 5I) 

^ (S61) 

if U{kx) for x S^) acts on 5”^ as a reflec¬ 

tion symmetry. By iterating Eqs. (S60) and (S61), any 
K-group in the present case reduces to that for a 1- 
dimensional nonsymmorphic insulator defined in the kx- 
direction. The Hamiltonian mapping (i=l—>-(7 = 2—>■ 
d = 3 discussed previously is based on the isomorphism 
(S60) and (S61). In Eq. (S60) or (S61), the first term 
in the right hand side represents a “weak” 
topological index of the left hand side, which is obtained 
by just neglecting the A^-dependence in the left hand 
side, but the second term gives the “strong” topological 
index. 


